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Abstract 

The Suzaku data of the highly variable magnetar IE 1547.0-5408, obtained during the 2009 
January activity, were reanalyzed. The 2.07 s pulsation, detected in the 15-40 keV HXD data, 
was found to exhibit phase modulation, which can be modeled by a sinusoid with a period of 
36 . 0 I 2 5 ks ar| d a n amplitude of 0.52 ±0.14 s. While the effect is also seen in the 10-14 keV 
XIS data, the modulation amplitude decreased towards lower energies, becoming consistent 
with 0 below 4 keV. After the case of 4U 0142+61, this makes the 2nd example of this kind 
of behavior detected from magnetars. The effect can be interpreted as a manifestation of 
torque-free precession of this magnetar, which is suggested to be prolately deformed under 
the presence of strong toroidal field of — 10 16 G. 

Key words: Stars:individual:1 E 1547.0-5408 — Stars:magnetars — Stars:magnetic field — 
Stars:neutron 


1 Introduction 

Magnetars (e.g., Thompson & Duncan 1995; Mereghetti 2008) 
are considered to be magnetically-powered isolated neutron 
stars (NSs) having ultra-strong dipole magnetic fields, Bd = 
1q14-i5 q as ,j er j ve( j f rom their pulse periods P and period 
derivatives P. These values of Bd exceed so-called critical 
field, 4.4 x 10 13 G, above which the Landau level separations of 
electrons exceed their rest-mass energy. Magnetars emit char¬ 


acteristic two-component X-ray spectra, consisting of the long- 
observed soft component and the mysterious hard component 
(Kuiper et al. 2006; den Hartog et al. 2006; Enoto et al. 2010b). 

Magnetars are thought to harbor even higher toroidal (or in¬ 
ternal) magnetic fields B t (e.g., Thompson & Duncan 1995; 
Harding & Lai 2006; Braithwaite 2009; DalLOsso et al. 2009). 
These fields may be generated and amplified in the progenitor’s 
core, during its final gravitational collapse involving strong dif- 
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ferential rotation (e.g., Takiwaki et al. 2009). Some models ex¬ 
plain the production of energetic particles in magnetars invok¬ 
ing toroidal magnetic fields and the associated twisted magnetic 
configuration (e.g., Beloborodov et al. 2013). 

Observationally, the presence of intense B t is supported by 
the recently discovered low-Bd magnetars, SGR 0418+5729 
(Rea et al. 2013), SGR1822-1606 (Rea et al. 2012), and 
3XMM J 185246.6+003317 (Rea et al. 2014); their burst activ¬ 
ity suggests that their internal fields are much stronger than their 
externally observable fields of Bd = 6 x 10 12 G, 2.7 x 10 13 G, 
and < 1.4 x 10 13 G, respectively. The possible proton cyclotron 
feature in SGR 0418+5729, which suggests a local surface field 
of > 2 x 10 14 G (Tiengo et al. 2013), provides another sup¬ 
port, because a part of such intense internal fields would emerge 
from the NS surface to form local multipoles. Thus, the strong 
B t may be regarded as one of the most important properties of 
magnetars. However, it remained difficult to more directly esti¬ 
mate B t , because toroidal fields are intrinsically confined inside 
the NS and invisible from outside. 

After pioneering works by, e.g., Carlini & Treves (1989) and 
Heyl & Hemquist (2002), the way around this difficulty was 
explored by Makisihma et al. (2014; hereafter MEA14) in the 
following way. When an NS has B t ~ 10 16 G, is it expected to 
be deformed into a prolate shape (Ioka 2001; Cutler 2002; Iokai 
& Sasaki 2004), to a degree of 

e={h- h)/h ~ 1 x 10- 4 (£ t /10 16 G ) 2 (1) 

where 1 3 is the moment of inertia around the NS’s symmetry 
axis x 3 , and Ii that around axes orthogonal to x 3 . Under the 
absence of external torque, such an NS has a constant rotation 
period P rot = 2nI 3 /L around x 3 , and a constant precession pe¬ 
riod Pprec = 2nh / L around the angular momentum vector L 
(with L = \L\) which is fixed to the inertial frame. These two 
periods, which were degenerate when spherical, thus become 
different just by e. If x 3 is tilted from L by a non-zero wob¬ 
bling angle a, the NS is said to exhibit free precession, which 
is the most basic behavior of a torque-free rigid body with an 
axial symmetry (Landau & Lifshitz 1976; Butikov 2006). This 
concept should be distinguished from the more familiar forced 
precession driven by external torque, as seen in spinning tops. 

Suppose that a precessing NS has a radiation pattern that 
breaks axial symmetry around * 3 . Then, as argued by MEA14, 
the observed pulses will be phase-modulated at the beat period 
between P rot and P prec , which is called slip period and is given 
as 

T = Pprec / C = (1/Prot - 1/Pprec ) _1 • (2) 

Using the Suzaku X-ray observatory (Mitsuda et al. 2007), 
MEA14 in fact found that the P = 8.687 s pulsation of the spec¬ 
tral hard component of the magnetar 4U 0142+61 is subject to 
phase modulation with a period of T = 55 ± 4 ks. When this 


effect is interpreted as a manifestation of free precession and T 
is identified with the associated slip period, equation ( 2 ) yields 
e = 1.6 x 10 -4 , and hence equation (1) implies B t ~ 10 16 G. 
This provides one of the first direct estimations of B t of a mag¬ 
netar. 

As seen above, the essential assumptions made by MEA14 
to interpret the observed phase modulation of 4U 0142+61 are; 
e ^ 0, a ^ 0, and the broken axial symmetry of the hard X- 
ray emission pattern around x 3 . If these assumptions generally 
hold in objects of this class, we should observe the same effects 
from other magnetars as well. If detected from a fair number 
of them, the effect will provide such important information as; 
how B t is distributed among them, whether Bd and B t are cor¬ 
related, and whether B t actually decreases (as postulated in the 
magnetar scenario itself; Thompson & Duncan 1995) towards 
objects with larger characteristic ages. The results will also pro¬ 
vide valuable clues to the origin of the hard X-ray component 
(Kuiper et al. 2006; den Hartog et al. 2006), to the nature of the 
intriguing spectral evolution of magnetars discovered by Enoto 
et al. ( 2010 b), and to the ultimate origin of these mysterious 
objects (e.g., Vink & Kuiper 2006; Nakano et al. 2015). Thus, 
examining other magnetars for the hard X-ray phase modulation 
is extremely important and urgent. 

In the present paper, we have chosen the highly variable 
magnetar IE 1547.0—5408 as our 2nd target. It is the fastest- 
rotating magnetar with P = 2.07 s, possibly associated with the 
supernova remnant G327.24—0.13 (Gelfand & Gaensler 2007). 
The observed value of P = 2.4 x 10 _11 s s _1 (Kuiper et al. 
2012), together with P, yields B d = 2.2 x 10 14 G and a very 
young characteristic age of r = P/2P = 1.4 kyr. On 2009 
lanuary 22, a strong bursting activity was detected from IE 
1547.0—5408 with Swift (Gronwall et al. 2009). This activ¬ 
ity was continuously covered by RXTE, INTEGRAL, and Swift 
observations (Kuiper et al. 2012). These monitoring studies re¬ 
vealed an increase in the total flux at the burst onset by more 
than an order of magnitude, followed by a rapid decline, and 
complex changes in the spin-down characteristics along the ac¬ 
tivity evolution. 

As already reported by Enoto et al. (2010a) (hereafter 
Paper I) and Enoto et al. (2012), this object was observed by 
Suzaku during its 2009 lanuary outburst, and yielded by far the 
strongest hard X-ray (> 10 keV) signals among some 10 mag¬ 
netars so far observed with Suzaku. Although this source was 
observed with Suzaku again on 2010 August, these follow-up 
data are not utilized here, since the source had faded by that 
time in hard X-rays to a level comparable to the HXD detection 
limit (Iwahashi et al. 2013). 
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2 Observation 

As described in Paper I and Enoto et al. (2012), the present 
Suzaku observation of IE 1547.0—5408 was performed from 
2009 January 28 21:34 (UT) to January 29 21:32 (UT), for 
a gross duration of 86 ks, with the target placed at the HXD 
nominal position. This was a Target-of-Opportunity observa¬ 
tion in response to the activity onset, and this epoch falls on 
“segment 3” (from 2009 January 22 to February 1) of Kuiper 
et al. (2012). Among the three operating cameras of the X-ray 
Imaging Spectrometer (XIS: Koyama et al. 2007), XIS0 was 
operated in the timing mode (P-sum mode), with a 7.8 ms time 
resolution. The Hard X-ray Detector (HXD: Takahashi et al. 
2007) was operated in the standard mode, yielding a time reso¬ 
lution of 61 ps. Further details of the observation are found in 
Paper I. 

Below, we mainly utilize the data from XIS0 and HXD-PIN 
processed as described in Paper I, in which each detected event 
is tagged with the energy and the arrival time. Those from XIS 1 
and XIS3 had too low a time resolution (2.0 s) to resolve the pul¬ 
sation, and signal statistics were rather low in HXD-GSO. The 
net exposures achieved with XIS0 was 42.6 ks, while that with 
the HXD was 33.5 ks. After subtracting the background and 
applying dead-time corrections, the object was detected with 
HXD-PIN at a 15-70 keV count rate of 0.174 ± 0.004 c s -1 
(Paper I), and a 1-10 keV rate of 5.5 ± 0.1 c s -1 with XIS0. 

As described in Paper I, 13 short-burst candidates were re¬ 
moved, and arrival times of the XIS0 and HXD-PIN events 
were photon-by-photon converted to those to be measured at 
the solar-system barycenter. The conversion utilized the source 
position of (a 2000 , 8 2000 ) = (15 h 50 m 54 B .ll,54 o 18'23."7), as 
well as the spacecraft orbital information. 

3 Data analysis and results 

3.1 Reconfirmation of Paper I 

Since the present paper re-analyzes the same data set of IE 
1547.0—5408 as used in Paper I, let us begin with reproduc¬ 
ing the timing results of that publication. While Paper I uti¬ 
lized the conventional chi-square method to search for the pul¬ 
sations, here we instead employ the Z 2 technique (Buccheri 
et al. 1983; Brazier 1994) which is free from the ambiguity 
of phase-bin number per pulse cycle and is generally less af¬ 
fected by the Poisson errors. For each trial pulse period P, this 
technique effectively calculates the folded pulse profile using a 
sufficiently large bin number (in practice not specified). Then, 
the Fourier power of the profile is summed up to a specified 
harmonic number n (n = 1,2,...), and normalized in a certain 
manner to reflect the data statistics. This process yields the sta¬ 
tistical quantity called Z 2 , which evaluates the significance of 
periodicity at P. If no intrinsic periodicity is present at P ex¬ 
cept Poissonian fluctuations, the Z 2 value calculated according 


to the normalization of Brazier (1994), which we adopt here, 
should obey a chi-square distribution with 2 n degrees of free¬ 
dom. The higher noise tolerance of this technique is realized by 
discarding the power in all harmonics higher than is specified. 

Figure la shows periodograms calculated in this way using 
the 2-10 keV XIS0 data and the 12-70 keV HXD-PIN data, 
both including backgrounds. The harmonic number was cho¬ 
sen to be n = 2, considering the double-peaked pulse profile of 
this source (e.g., Paper I). These energy bands, as well as the 
trial period range (2.0713-2.0730 s), are the same as those in 
Paper I, while the employed period search step, 1 x 10~ 6 s, is 4 
times finer. The results reconfirm essential features of figure 3 
of Paper I, including the best-estimate period of 

P 0s = 2,072.135 ±0.005 ms (3) 

referring to the XIS result, and the side lobe patterns. The HXD 
periodogram peak, 

P oh = 2,072.147 ±0.005 ms, (4) 

is consistent with equation (3) within errors, although it is 
slightly displaced to longer periods. This behavior is also seen 
in figure 3 of Paper I. The errors in Po s and P 0 h refer to the 
range where pulse-phase differences amount to ±10% of one 
pulse cycle when accumulated over the gross duration of 86 ks. 
These errors approximately agree with widths of the peaks in 
the periodograms. 

At P oh , the 12-70 keV HXD periodogram in figure la 
reaches the peak height of X = 28.21, where Zf is abbreviated 
to X. In terms of chi-square distribution with 2n = 4 degrees 
of freedom, the chance probability of finding X values higher 
than this is found as p£ = 1.1 x 10 -4 ; here, the subscript spec¬ 
ifies the harmonic number n, and the superscript “h“means the 
HXD. This agrees with the probablity of ~ 1 x 10 -4 obtained 
in Paper I in 10-70 keV with the chi-square technique. 

In figure la, the XIS and HXD periodograms both show a 
number of side lobes. Such a structure at a period P' often 
arises due to beat between the pulsation at Po (representing Po s 
or P oh ) and a certain longer period T. Based on this assumption, 
we calculated the required period as T = |(l/Po — 1 /P') - 1 1, 
and show it in green at the top of figure 1. In addition, the 
required values of T (in units of ks) of major XIS0 and HXD 
side lobes are given explicitly in the figure. Thus, many of the 
side lobes (particularly of XIS0) in figure la can be explained 
as a beat of Po with T = 5 — 6 ks or T = 10 — 12 ks, which 
agree with the Suzaku's orbital period P or b = 5.6 ks or twice 
this value, respectively. However, the HXD-PIN periodogram 
shows several unexplained peaks; at T = 24.8 ks and 38.7 ks on 
the left side of Po, as well as at T = 51.8 ks, 21.1 ks, and 18.1 
ks on the right side. These structures, also seen in figure 3 of 
Paper I, suggest that the hard X-ray pulses could be affected by 
some long periodicity in the range of T = 18 — 52 ks (5-14 hr). 
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3.2 Detailed analysis of the HXD data 

Below, we analyze the HXD-PIN data in further details, select¬ 
ing an energy range of 15-40 keV. The lower energy bound is 
thus raised from 12 keV used in figure la to 15 keV, to avoid 
increased soft thermal noise. Likewise, the upper energy bound 
is lowered from 70 keV to 40 keV, to reduce the contribution 
from the non X-ray background. 

3.2.1 Evolution of hard X-ray pulse profiles 

The Z\ periodogram calculated from the 15-40 keV HXD-PIN 
data is shown in figure lb in black. Except some small differ¬ 
ences, it is similar to the 12-70 keV results shown in brown in 
panel (a). The peak value, also obtained at the period of equa¬ 
tion (4), becomes 

A' 0 = (X ) A =o = 22.06 , (5) 

which is somewhat smaller than that in 12-70 keV. 

In the 15—40 keV periodogram, the period range of T = 
18 — 52 ks, where the unexplained side lobes were noticed in 
figure la, is approximately preserved (although detailed of the 
side lobes changed). Since this period range does not contain 
particular artificial periodicities of the observation such as data 
gaps or background variations, some properties (e.g., amplitude, 
phase, or pulse shapes) of the 2.07 s hard X-ray pulsation could 
be subject to slow variations with periods in this range. In fact, 
Kuiper et al. (2012) found that the pulse shapes were change¬ 
able during their “segment 3” which covers the present Suzaku 
observation, that a new intermediate peak in the > 4 keV pulse 
profile emerged, and that the pulsed hard X-ray flux increased 
meantime although the total 4-27 keV flux was rapidly declin¬ 
ing. 

To examine the evolution of the hard X-ray pulse proper¬ 
ties, we divided the whole observation into consecutive time 
segments of approximately equal lengths, denoted as j = 
0,1,..., J — 1 (j = 0 at the beginning and j = J at the end), and 
sorted the entire 15—40 keV HXD-PIN events (including back¬ 
grounds) into a two-dimensional array C(j,k). Here, the col¬ 
umn k = 0,1.., K — 1 represents the pulse phase (i.e., the photon 
arrival times modulo Poh). Thus, a set {C(j,k)\k = 0,..,K — 1} 
gives the pulse profile in the j- th time segment. We employed 
J — 6 and K = 12, and show in figure 2a the obtained time- 
divided six pulse profiles as a function of k. There, running av¬ 
erages were applied to 3 consecutive data points (k — l,k,k + l). 
Thus, the pulse amplitude changes with time, and the pulse 
shape is also variable; it is relatively flat-topped with a single 
deep minimum, with an occasional shallower secondary mini¬ 
mum at about half a cycle offset. Such a short-term variation 
in the pulse shape was observed from 4U 0142+61 (Enoto et 
al. 2011), as well as from IE 1547—5408 itself (Kuiper et al. 
2012 ). 

More importantly, thanks to the high sensitivity of the 


Suzaku HXD on relatively short time scales (Kokubun et al. 
2007), the pulse phase is seen to vary clearly in figure 2a, by 
up to about half a pulse cycle. On the other hand, we do not 
notice in figure 2a particular secular changes in the pulse pro¬ 
file, pulsed flux, or pulse phase. This does not contradict to the 
discovery by Kuiper et al. (2012) of the emerging new compo¬ 
nent in the hard X-ray pulse profile, because the growth of that 
component is estimated to be at most ~ 10% across the present 
Suzaku pointing, as judged from the history of the hard X-ray 
pulsed flux (figure 9 of Kuiper et al. 2012). 

3.2.2 Periodic phase modulation 

We next need to know whether the pulse-phase fluctuations seen 
in figure 2a are periodic or not. This requires a finer time 
division by increasing J, but individual pulse profiles would 
then lose statistics. Therefore, we chose another way; the 
data accumulation was repeated with J = 20 and K = 5, but 
the results are presented in figure 3a in the form of five light 
curves, {C(j, k)\j = 0,.., J = 19} using j as abscissa. For pre¬ 
sentation, 3-point running averages were again applied over 
(j — 1 ,j,j + 1). Thus, the background-inclusive count rates in 
the five pulse phases mostly varied rather periodically, by 2 to 3 
cycles across the 86 ks gross pointing, although the amplitude 
is comparable to the statistical errors. 

We calculated discrete Fourier power spectra of the five light 
curves in figure 3a, before taking the running averages, and 
summed them to obtain a single power spectrum which is shown 
by a red solid line in figure 3c. It indeed reconfirms the sug¬ 
gested power enhancement, at wave numbers of q = 2 and q = 3, 
with the corresponding period range of 

86/3 = 29 ks to 86/2 = 43 ks . (6) 

The weighted mean of q = 2 and q = 3 over the power spectrum 
in figure 3c becomes 

q — 2.45, or T = 35.Iks, (7) 

which is fully consistent with the implication of the side lobes 
in figure 1. 

As a “control” study, figure 3b presents another set of light 
curves, calculated in the same way but purposely changing the 
folding period from Po to a dummy value of 2.0 s. Its behav¬ 
ior is clearly more random, without particularly favored peri¬ 
ods. Indeed, these 5 control light curves yielded the blue power 
spectrum in figure 3c, which is consistent with white noise. 

To evaluate the significance of the q = 2 and q — 3 enhance¬ 
ments in the red power spectrum of figure 3c, we repeated the 
same control study 20 times, by changing the dummy period 
around 2.0 s. The derived ensemble-averaged power spectrum, 
shown in figure 3c by a dashed black line, is approximately 
white (with hints of red noise at q = 1 and background variation 
at q = 4 corresponding to 5P or b). Around this mean, individual 
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power spectra exhibited a lcr scatter by ~ ±0.5. Therefore, the 
q = 2 and q = 3 excess (above 1.0) in the red power spectrum is 
estimated to be 2.4 and 1.6 sigma effects, and their sum implies 
a ~ 3 sigma excess. 

Since the five light curves in figure 3a do not vary in phase, 
the suggested long periodicity is unlikely to result from pulse- 
amplitude variations, but considered to reflect the pulse-phase 
jittering found in figure 2a. Thus, the hard X-ray pulsation 
is inferred to exhibit periodic phase modulation. To reinforce 
this inference, we accumulated the pulse profiles in 6 separate 
phases of an assumed periodicity of T = 36 ks, of which the 
choice is justified later. The result is presented in figure 2b; in 
this case, the row j of the array C(j, k) denotes the pulse phase, 
i.e., the arrival times modulo T, rather than the arrival times 
themselves. It visualize a periodic and approximately sinusoidal 
pulse-phase modulation, with an amplitude of A k ~ ±3.5 bins 
(±0.3 pulse cycle or ±0.6 s) which is consistent with the infer¬ 
ence from figure 2a. 

3.2.3 Demodulation analysis 

The above semi-quantitative analysis justifies us to apply the 
“demodulation analysis”, developed by MEA14, to the HXD 
data, to more rigorously quantify the pulse-phase modula¬ 
tion. Specifically, we consider that the hard X-ray pulsa¬ 
tion of IE 1547.0—5408 is somehow phase modulated, like in 
4U 0142+61, at a period T in the range of equation (6). We 
hence hence assume that the peak timing t of each 2.07 s pulse 
is displaced sinusoidally by 

At = Asin(2nt/T — (fi) : (8) 

where A and 0 are the amplitude and initial phase of the as¬ 
sumed modulation, respectively. 

As formulated in MEA14, such effects can be removed 
by shifting the arrival times of individual HXD photons back 
by —At. This procedure was previously employed, e.g., by 
Koyama et al. (1989) to search for orbital motions in the anoma¬ 
lous X-ray pulsar IE 2259+586, and by Terada et al. (2008) 
to remove the orbital effects in the fastest-spinning accreting 
white dwarf binary, AE Aqr, and successfully detect its hard 
X-ray pulsations. We hence applied these time displacements 
to the individual HXD photons, and re-calculated the Z | peri- 
odograms to see whether the pulse significance changes. The 
pulse period was scanned over a range of P = 2.072130 — 
2.072152 (with a step 1 /rsec), i.e., the joint error range of equa¬ 
tion (3) and equation (4). The triplet (T, A, 0) describing the 
phase modulation was scanned over A = 0 — 0.80 sec (with a 
0.05 sec step), 0 = 0- 360° (3° - 10° step), and T= 10-80 
ks (with a step of 0.2 to 0.5 ks). The harmonic number n = 2 is 
retained. 

The results of this demodulation analysis are presented in 
figure 4, just in the same format as figure 2 of MEA14. There, 


panel (a) shows a two-dimensional color map of the value of 
X = Zl (the maximum as P is varied) obtained on each grid 
point (A, 0), for a particular case of T = 36.0 ks. Its projections 
on the 0 and A axes are presented in panels (b) and (c), respec¬ 
tively, where multiple traces indicate cross sections of panel (a) 
at various positions. While the value of Xo in equation (5) is re¬ 
produced at A = 0 in panel (c), a clear peak is seen at A ~ 0.52 
s and 0 ~ 270°. This value of A agrees well with what can be 
read from figure 2b. Furthermore, the value of 0, together with 
equation (8), predicts At = +A (the maximum delay) at about 
the beginning ( t ~ 0) of the T = 36 ks phase, again in agreement 
with figure 2b. Although figure 2a apparently suggests a rather 
long periodicity (e.g., T > 100 ks), this is probably a beat be¬ 
tween the T = 36 ks modulation and the Nyquist period therein 
(1 x 86/6 = 29 ks), which should appear at 140 ks. 

In panels (a)-(c) of figure 4, the peak significance reaches 
Xp k = 45.60; here, X p k(T) denotes the peak value of X found 
on a (A,0) plane for a given value of T. Panel (d) of figure 4 
shows how this X p k (T) depends on the assumed modulation 
period T. Over the entire (P,T, A,0) space surveyed, the high¬ 
est value of X, to be called the grand-maximum and denoted 
Xmax, is thus found as 

X max = X p k(36.0ks) = 45.60 , (9) 

for a set of modulation parameters given as 

A = 0.52 ± 0.14s, 0 = 276° ±15°, T = 36.0±H ks . (10) 

Here, the errors are represented in the same way as in MEA14, 
i.e., by the standard deviations of Gaussians fitted to the dis¬ 
tributions (above uniform backgrounds) in figure 4. This value 
of T, already employed to produce figure 2b, agrees well with 
equation (7). 

As shown in figure lb in red, the demodulation with equation 
(10) has changed the HXD-PIN periodogram significantly. The 
peak at Poh has become much higher, just from Xo to X max . 
Furthermore, the P or b-unrelated side lobes have diminished; 
powers therein are thought to have been restored to the main 
peak. 

In figure 2c, we again calculated the 15^-0 keV folded pulse 
profiles in the same 6 time segments as in figure 2a, but af¬ 
ter applying the time corrections by equation (8) to individual 
photons using the parameters in equation (10). Thus, the con¬ 
secutive 6 profiles have indeed lined up in their phases. A small 
subpeak, just prior to the minimum, is seen in at least 4 out of 
the 6 profiles. These results reinforce the validity of the demod¬ 
ulation analysis. Figure 2d reproduces, and directly compares, 
the time-averaged 15-40 keV pulse profiles before and after the 
demodulation. Thus, the pulse amplitude has been nearly dou¬ 
bled. 



6 


Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 


3.2.4 Significance of the phase modulation 

Although the phase demodulation procedure largely increased 
the HXD pulse amplitude (and hence the pulse significance), we 
need to examine whether this increase is statistically significant, 
or merely due to a chance superposition of Poisson noise of the 
data on the already existing pulse periodicity at Poh- 

The nominal chance probability of finding values of X > 
A'max of equation (9) is extremely small, = 3.0 x 10~ 9 . 
However, we cannot use this face value to evaluate the statis¬ 
tical significance of the detected effect, for the following two 
reasons. 

1. This must be multiplied by the total number of inde¬ 
pendent trials involved in the computation of figure 4, but 
this quantity is not easy to estimate. The total number of 
parameter search steps, actually employed, would consider¬ 
ably overestimate this number, because the steps were cho¬ 
sen, rather arbitrarily, to be fine enough not to miss the X 
peaks. 

2. When evaluating the significance of X max , we need to con¬ 
sider the fact that the HXD signals were already pulsing, be¬ 
fore the demodulation procedure, at Poh (Paper I) with a sig¬ 
nificance of Xo n. The variations in the pulse amplitude, 
profiles, and phase (figure 2) would further contribute to X. 

These two problems can be combined into a single question: 
when we scan over P, A, <j> and T, like in figure 4, what is 
the overall chance probability for the X statistics to increase by 
chance from Xo to X max or higher when the source is pulsating 
with variable parameters but without any excess variation at T 
above the white noise? 

To answer the above inquiry, we conducted 1,000 simulation 
runs as detailed in Appendix, using the actual 15—40 keV data 
themselves, but purposely randomizing, each time, the phase of 
the assumed 36 ks modulation: any white noise would not be 
affected by such randomization. Figure 5a shows a distribution 
of the grand-maximum value X max obtained in this simulation 
study. It is considered to represent the probability distribution 
of A max , to be found in the same demodulation analysis of the 
HXD data over the same ( P, T, A, rj>) space as figure 4, but un¬ 
der the absence of power enhancement (above a white noise) 
at any T examined. There, values exceeding equation (9) were 
actually found in three out of the 1,000 runs, namely, with a 
probability of 0.3%. 

Figure 5b gives the upper integral probability associated 
with figure 5a. In the range of X max > 35, it is approximated 
by an exponential as 


/(> X m ax) = /oexp 


X max - 40 
9o 


(ID 


with fo = 0.038 ± 0.005 and go = 2.2 ± 0.2, where X max = 
40.0 was chosen as a reference pivot point. The value of go 
is consistent with the exponential factor, 2 . 0 , that is expected 


for a chi-square distribution of any degree of freedom. Thus, a 
value of X max which is equal to or larger than equation (9) will 
appear with a chance probability of 

/(45.60) = (3±1) x 10 -3 , (12) 

in a good agreement with the implication of panel (a). As no¬ 
ticed in Appendix, this estimate is considered rather conserva¬ 
tive. 

Considering that MEA14 used n = 4, it may be interest¬ 
ing to examine what happens if different values of n are em¬ 
ployed. We hence repeated the same analysis with n = 1, 3, 
and 4, to find X max = 25.64, 47.06, and 49.92, respectively. 
The optimum parameters were not much different from those 
in equation (10). Although the larger harmonic numbers some¬ 
what increased X max from equation (9), the effect is not sig¬ 
nificant considering the increased degree of freedom. In fact, 
the implied chance occurrence probability, = 2.8 x 10 -6 , 
V\ = 1.8 x 10~ 8 , and V 4 = 4.2 x 10 -8 , are all much higher 
than V 2 obtained above, in agreement with the double-peaked 
pulse profiles. Therefore, we retain n = 2 below. 

In short, we can exclude, at > 99.6% confidence, the possi¬ 
bility that Xmax of equation (9) in figure 4 arised via chance 
superposition of the photon count fluctuations on top of the 
hard X-ray pulsation of which the parameters vary with white 
noise. This confidence level is more stringent than the case of 
4U 0142+61, where it was > 99% (MEA14). 

3.3 Analysis of the XIS data 

The two-component spectral model, fitted in Paper I to the 
XIS+HXD spectrum of IE 1547.0—5408, indicates that the soft 
thermal component and the hard power-law component cross 
over at ~ 5 keV. Therefore, the hardest end of the XIS energy 
band, with its P-sum time resolution, is expected to provide in¬ 
dependent information on the hard component, and allows us 
to crosscheck the HXD-PIN results. (This was not feasible with 
4U 0142+61, since the XIS was not operated in the P-sum mode 
in that observation.) To limit the soft-component contribution to 
< 10%, we choose an energy range of 10-14 keV, and perform 
the same analysis as above. Although this energy range, which 
is in fact covered by the XIS, is seldom utilized in spectral stud¬ 
ies mainly because of the calibration uncertainty therein, such 
problems do not affect timing studies. 

Figure 6 show the results of this demodulation analysis on 
the 10-14 keV XIS0 data, presented in the same way as figure 4. 
Panels (a), (b), and (c) all refer to the condition of T = 36.0 ks 
after equation (10). On the (A, <j>) plane, a clear peak with 

X pk (36ks) = 27.64 (13) 

is seen at a similar position to figure 4a. More specifically, A 
and <j> can be constrained as 
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A = 0.46 ± 0.11s, 0 = 260° ±17° . (14) 

These values agree, within errors, with the corresponding re¬ 
sults in equation (10). If the modulation parameters are fixed to 
the best-fit values of equation (10) as determined with the HXD, 
we instead find 

X pk (36ks) = 17.08 . (15) 

The demodulation with equation (14) has effects on the 10- 
14 keV periodograms as shown in figure 7a. While the value 
of X ~ 5 found before the demodulation (black) is within ran¬ 
dom fluctuations, a clear peak, of which the height is given by 
equation (13), appeared at a period that is in full agreement with 
equation (3) and equation (4). 

Figure 7b shows 10-14 keV XIS0 pulse profiles folded at 
Pos, obtained using the demodulation parameters of equation 
(14) and presented in the same colors as in figure 2d. Again, the 
demodulation procedure significantly increased the pulse frac¬ 
tion. The red pulse profile is similar to that in figure 2d, al¬ 
though the relative depths of the two pulse minima are reversed. 
For reference, 1-10 keV XIS0 pulse profile is also shown. 

The T-dependence of X pk is shown in figure 6 d. Thus, in 
addition to a small hump at T — 36 — 38 ks, we observe many 
other peaks, some being much higher. Among them, the sharp 
feature at T ~ 11 ks is likely to be caused by background vari¬ 
ations, because it is close to 2P or b. Besides, the plot reveals 
a prominent peak at T = 44 ks, which is probably due to a 
chance superposition of Poisson fluctuations on top of the de¬ 
modulated pulsation. As a result, we cannot claim the pres¬ 
ence of the T = 36 ks periodicity based on the XIS data alone. 
Nevertheless, the XIS results significantly reinforce the HXD- 
detected phase modulation: since the pulsation is insignificant 
in the raw periodogram (figure 7a), the probability of finding a 
peak of Xpk = 17.08 of equation (15) at the same ( P,T,A ,</>) 
point as specified by the HXD result, can be directly evaluated 
with a simple chi-square distribution to be "P| = 1.9 x 10~ 3 . If 
we instead use equation (13), the probability further decreases 
to VI = 1.5 x 10 -5 . 

3.4 Energy dependence of the modulation amplitude 

Now that the 36 ks phase modulation in the hard X-ray pulses, 
first detected with the HXD, has been reinforced by the 10-14 
keV XIS data, our next task is to examine how this phenomenon 
behaves towards lower energies. In the same way as above, we 
hence analyzed the XIS0 data in three lower energies; 7-10 keV 
where the spectral hard component is still dominant, 4-7 keV 
where they compete, and 1=1 keV dominated by the soft com¬ 
ponent. Figure 8 shows the results of this soft X-ray analysis, 
namely, scans in A in which T — 36.0 ks is fixed while <j> is al¬ 
lowed to take any value. Thus, the 36 ks modulation amplitude 
deary decreases towards lower energies, and the effect almost 


disappears in the 1-4 keV range. 

Taking the 7-10 keV case as an example, we again calcu¬ 
lated the 12-bin folded pulse profiles for 6 phases of T = 36 ks, 
in a similar manner to figure 2b. The result, shown in figure 8d, 
again visualizes the sinusoidal phase modulation, clearly with a 
smaller amplitude than in figure 2b. The amplitude can be read 
as A k/K ~ ±2/12 = ±0.17 pulse phase, or A = 0.17Poh = 0.35 
s, in a good agreement with figure 8a. In this case, the pulse 
peak is more clearly traced than the valley. Furthermore, the 
modulation phase is somewhat shifted from equation (10) to 
4> ~ 360°. 

Combining figure 8 with equation (10) and equation (14), 
the energy dependence of A is visualized in figure 9. There, the 
dashed curve represents the fractional contribution of the hard 
component at that energy, namely 

A = A 0 xF h /(F a + F h ) (16) 

where Ao = 0.5 s is a fiducial amplitude at the hardest energy 
limit, while F s and Fh are the intensities of the soft and hard 
components, respectively, according to the spectral decompo¬ 
sition (Model B) in Paper I. Interesting enough, this empirical 
quantity behaves very similarly to the observed energy depen¬ 
dence of A, suggesting that the phase modulation is limited to 
the hard X-ray component. 

Given figure 9, we attempted to divide the 15—40 keV HXD- 
PIN band into finer ranges, or to extend the analysis to higher 
energies beyond 40 keV up to ~ 70 keV. However, due to lim¬ 
ited statistics, the data provided no significant information on 
the possible energy dependence of the HXD phase modulation. 
Likewise, the HXD-PIN data below 15 keV (down to 10 keV) 
had so poor statistics that they did not allow us to cross con¬ 
firm the 10-14 keV XIS results. All what can be said is that 
the pulse-phase modulation in the HXD-PIN signal is consis¬ 
tent with being energy independent over the 10-70 keV range 
within rather large uncertainties. 

3.5 Analysis of the HXD-GSO data 

As described in Paper I, the signal from IE 1547.0—5408 was 
marginally detected in the 60-100 keV range with HXD-GSO. 
Therefore, we analyzed the HXD-GSO data in various energy 
ranges from 50 to 200 keV, but the data gave significant evi¬ 
dence for neither the pulsation, nor its phase modulation. This 
is not surprising, because HXD-GSO has a considerably higher 
background than HXD-PIN (Kokubun et al. 2007). 

Even though the HXD-GSO data provide no information 
on the source pulsation, they allow us to perform an impor¬ 
tant control study. We hence purposely chose a very high en¬ 
ergy range, 300-307 keV, which is dominated by background 
and has nearly the same counting rate, 0.46 ct s -1 , as the 
background-inclusive 15-40 keV HXD-PIN rate (0.44 ct s -1 ; 
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figure 2). The green trace in figure 4d shows the derived X p k 
behavior in this high energy band. Thus, the baseline value of 
A'pk ~ 10 is much lower than that in 15^-0 keV, A pk ~ 22, 
due to the absence of the pulsation. In addition, no peaks 
are observed around T ~ 36 ks. Since these HXD-GSO data 
were acquired by the same instrument as the HXD-PIN data, 
at completely the same time, and under the same background 
conditions, the result argues against instrumental origins of the 
T = 36 ks periodicity in the 15-40 keV pulse phases. 

4 Discussion 

Through an extensive reanalysis of the Suzaku data of IE 
1547.0—5408 acquired in the 2009 outburst, we obtained ev¬ 
idence that the 2.07 s hard X-ray pulsation detected with the 
HXD is subject to phase modulation with a period of T = 36 
ks (figure 4). Following the first detection from 4U 0142+61 
(Po = 8.69 s, T = 55 ks, and A ~ 0.75 s; MEA14), this makes 
IE 1547.0—5408 a 2nd magnetar exhibiting this kind of evi¬ 
dence. 

4.1 Trivial possibilities 

4.1.1 Instrumental origin 

The detected effect should be first examined for instrumen¬ 
tal or observational artifacts. In section 3.2.4, we excluded 
the statistical fluctuation origin at a 99.6% confidence level. 
Furthermore, the HXD-GSO data analysis (section 3.5) indi¬ 
cates that the HXD instrument does not suffer unexpected ar¬ 
tifacts (e.g., background variations) that can explain away the 
A' p k maximum at T = 36 ks. The value of T = 36.0 ks 
= 6.4P or b is not in harmonic resonance with the Suzaku’s or¬ 
bital period, either. 

In this respect, of particular importance is the cross confir¬ 
mation of the phase modulation with the 10-14 keV XIS data 
(figure 6). This further rules out the possibility of the effect 
being due to some artifact specific to the HXD. Furthermore, 
the clear decrease of A towards lower energies (figure 9) ex¬ 
cludes any instrumental or artificial origin specific to the overall 
Suzaku mission (e.g., unexpectedly large drifts in the spacecraft 
clock), because such effects should be independent of the pho¬ 
ton energy. We are hence left with celestial interpretations. 

4.1.2 Period derivatives 

Based on figure 2b, figure 3a, and figure 3c, we have so far re¬ 
garded the phase modulation as periodic, and performed the de¬ 
modulation analysis to arrive at the modulation period of equa¬ 
tion (10). However, the total data span, 86 ks, contains only 2.4 
cycles of the 36 ks oscillation. Hence, we could be mistaking 
some other effects for the periodic phase changes. 

One possibility is the period derivative, because magnetars 


generally spin down, presumably emitting magnetic dipole ra¬ 
diation. This will cause parabolic changes in the pulse phase, 
which could mimic a periodic phase modulation. From Entry 
No.3 of Table 4 in Kuiper et al. (2012), the average spin-down 
rate of IE 1547.0—5408 is considered to be P ~ 2 x 10~ 11 s 
s ' 1 around the present observation. Although this would cause 
a phase difference of only ^ 0.1 pulse cycle across the 86 ks 
length of the present data, there still remains a possibility that 
the instantaneous P was considerably higher. 

With the above in mind, the Z\ analysis was repeated on the 
15-40 keV HXD-PIN data, incorporating this time the P term. 
We defined the period Po as the value at the middle point of the 
observation, and scanned it over the same joint error range as 
defined by equation (3) and equation (4). The analysis was first 
conducted assuming no phase modulation (A = 0). However, 
the pulse significance did not increase significantly; instead, the 
X value decreased from equation (5) as |P| is increased. Thus, 
a loose limit was obtained as —1.5 x 10 -10 < P < 3.2 x 10~ 10 
s s -1 . We next fixed T = 36 ks, and examined the parameter 
region around equation (10). However, the pulse significance 
did not increase, either. The obtained limit is again uninterest¬ 
ing, — 1.1 x 10~ 10 < P < 0.8 x 10 -10 s s -1 . In short, period 
derivatives cannot explain the observed pulse-phase modulation 
effects. 

4.1.3 Non-periodic origin 

Another possibility is that the hard X-ray pulse profile (includ¬ 
ing its phase) is fluctuating randomly for some unspecified rea¬ 
sons, and the variability power happened to be high at T ~ 36 ks. 
To help examination of this case, the wave number q = 86 ks /T 
is drawn in blue at the top of figure 4d. Thus, the significant ex¬ 
cess in A'pk = Z\ above equation (5), e.g., A pk > 30, is confined 
to a rather narrow range of the wave number as q = 2.4 ± 0.5. in 
a good agreement with figure 3c and equation (7). Furthermore, 
we have already confirmed in section 3.2.2 that the power en¬ 
hancements at q = 2 and q = 3 in figure 3c are 3 sigma effects. 
These facts altogether mean that the phenomenon should be re¬ 
garded as periodic rather than random. 

There is additional evidence arguing against random phase 
fluctuations, at least 10 < T < 30 ks is concerned. In this range, 
the 15—40 keV A p k values in figure 4d are comparable to Xq, 
and are all found at A < 0.15 s, implying that the pulse signifi¬ 
cance diminishes by demodulation with T in this range. Over a 
longer period range of 30 < T < 50 ks, A pk in figure 4d is ob¬ 
tained at A = 0.5 — 0.6 s, as a simple continuation from equation 
(10). At T > 50 ks, the values of A become rather poorly con¬ 
strained because of the insufficient data length compared to T. 
Nevertheless, the pulse-phase modulation is not likely to con¬ 
tain significant power at these long periods, either, because the 
power spectrum in figure 3c shown in red is not enhanced at 
q — 1 above the Poissonian contribution. Therefore, the pulse- 
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phase modulation is unlikely to have “pink noise” nature includ¬ 
ing secular changes, which would usually contain significant 
power at q = 1 

To assess the above issue yet from another aspect, we ap¬ 
plied the demodulation analysis separately to the first and sec¬ 
ond halves of the 15-40 keV HXD data. The procedure is the 
same as in subsection 3.2.3, except that T = 36.0 ks was kept 
fixed because of the shorter data length. Then, the first half 
gave the significance peak with X max = A' p k(36) = 32.4 at 
A = 0.58 ±0.18 s and <f> = 272° ± 18°, while the second half 
yielded X max = 22.7 at A = 0.42 ± 0.17 s and 0 = 293° ± 30°. 
In terms of the parameter (A,4>), the two halves agree within 
errors with each other, and with equation (10). Thus, the 36 ks 
pulse-phase modulation is considered to have persisted through 
the present observation. 

Finally, figure 10 shows the 15-40 keV pulse profiles in the 
1 st and 2 nd halves of the data, demodulated using the common 
set of parameters given by equation (10). In fact, we discarded 
the initial 14 ks of the 1 st half and the final 14 ks of the 2 nd 
half, to make them both represent exactly the time period of 36 
ks. Within the photon count fluctuations, the two profile agree 
with each other, with \ 2 = 0-78 for v = 12 (calculated before 
applying the running average). Thus, the pulse profiles from 
the 1 st and 2 nd modulation cycles, when demodulated using the 
common set of parameters, become consistent with each other. 

In conclusion, we can regard the phase variations in the hard 
X-ray pulses as periodic, in spite of the rather limited data span. 

4.2 Interpretation as free precession 

4.2.1 Free precession 

Let us hereafter limit our interpretation to periodic celestial ef¬ 
fects. The simplest of them is binary motion of the magnetar: 
indeed, A = 0.5 s and T = 36 ks can be explained if the puta¬ 
tive companion star has a mass of the order of 0.1 solar mass 
(depending on the inclination). However, even putting aside the 
astrophysical reality of such a binary, the absence of the phase 
modulation in soft X-rays (figure 9) clearly rules out this inter¬ 
pretation, just as in 4U 0142+61. 

From the above arguments, free precession (Carlini & Treves 
1989; Heyl & Hernquist 2002; MEA14) described in section 1 
remains the most promising possibility. In this scenario, we 
identify x 3 with the dipole magnetic axis, assume that it is tilted 
from L (i.e., a ^ 0), and consider that the NS is deformed to 
have a non-zero value of e = A I/I. Then, P TO t and P prec , 
which were degenerate when e = 0 , split into slightly different 
two values according to equation ( 1 ). Of them, P prec around 
L is observed as the pulse period of equation (3) and equation 
(4), i.e., P P rec = Po, while P TOt around x 3 cannot be detected 
when the emission pattern is symmetric around x 3 even if t ^ 0 
and a^0 both hold. This case is thought to apply to the soft 


component, because it is presumably thermal emission from the 
NS surface, and should be symmetric around the magnetic axis 
along which the heat will flow. 

If, in contrast, the emission pattern breaks symmetry around 
* 3 , the observed signals would depend not only on the direction 
of £3 as seen from us, but also on the NS’s rotation angle around 
£3 relative to the L-x 3 plane. Since this relative angle changes 
with the slip period T of equation (2), the pulse phase becomes 
modulated at T (MEA14). With a natural identification of T 
with the modulation period of 36 ks, this case can explain the 
behavior of the hard X-ray component, and hence equation (2) 
yields 

N = |A///| = P 0 /T = (0.58tS$) x 10~ 4 . (17) 

4.2.2 Double periodicity 

If the above interpretation is correct, we should be able to ob¬ 
serve periodicity also at P rot = P pr ec/(1 + e). From equation 
(17), it should appear as a prograde period at 

P+ t = 2,072.027 ±0.014 ms (18) 

if the deformation is prolate (e > 0 ), and a retrograde period at 

P~t =2,072.267 ±0.014 ms (19) 

if oblate (e < 0 ), although they could appear simultaneously. 
Of them, P rpt is already visible in figure lb. Nevertheless, the 
issue deserves further examination, because the periodicity at 
Prot would also be subject to the phase modulation at T, and 
the significance would change by the demodulation procedure. 

Figure 11 shows all the periodograms obtained while com¬ 
puting figure 4, superposed onto a single plot. In addition to 
the main peak coincident with equation (4), the results thus re¬ 
veal clear enhancements at ~ 2,072.03 ms and ~ 2,072.25 ms. 
Since they can be readily identified with P+ t and P“ t , respec¬ 
tively, our expectation has been confirmed. Furthermore, the 
P+ t peak is somewhat more significant, in agreement with our 
assumption that the object has a prolate form. To more quan¬ 
titatively asses the significance of these two enhancements and 
their differences remains our future study. 

4.3 Comparison with 4U 0142+41 

In terms of free precession of an axisymmetric NS, let us com¬ 
pare IE 1547.0-5408 and 4U 0142+41 (MEA14). The most 
important feature common to them is the similar degree of de¬ 
formation, t=1.6x 10 -4 in 4U 0142+41 and t = 0.6 x 10 -4 
IE 1547.0—5408 [equation (17)]. This suggests that the de¬ 
formation in the two objects has a common origin. Another 
resemblance between them is the decrease of the pulse-phase 
modulation amplitude towards lower energies; at least in IE 
1547.0—5408, the effect is consistent with being absent in the 
soft spectral component. As already described in section 4.2.1, 
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these properties can be interpreted by considering that the soft 
component is emitted symmetrically around the stellar magnetic 
axis, while the hard component breaks that symmetry. This dif¬ 
ference suggests that the two spectral components are distinct 
not only in their spectral shapes, but also in their emission pat¬ 
terns and/or emission regions. 

While the two sources thus show similar behavior, they also 
exhibit the following three differences, which must be explained 
in order to further solidify the interpretation. 

1. While we measured A/Po =0.08 ±0.04 in4U 0142+61, the 
ratio is much higher, A/Po = 0.25 ±0.07 (approximately 
±l/4pulse cycle) in IE 1547.0—4508 in energies above 10 
keV (figure 9). 

2. While n = 4 was needed in 4U 0142+61, n = 2 is optimum 
in IE 1547.0-45408 (section 3.2.4). 

3. The hard X-ray pulsation of 4U 0142+61 in the 2009 
Suzaku observation was detected only after the demodula¬ 
tion, whereas that of IE 1547.0—5408 is significant even be¬ 
fore such corrections (Paper I). 

Let us consider item No.l first. In 4U 0142+61, the sinu¬ 
soidal hard X-ray modulation with A/Po = 0.08 was attributed 
to positional (or directional) deviations of the X-ray emission 
region (or direction) from the magnetic poles (or from £ 3 ). 
There, the emission pattern was implicitly and without particu¬ 
lar reason assumed to be pencil-beam like. In contrast, the case 
of IE 1547.0—5408 can be better explained by assuming that 
the emission reaches us (for simplicity, from a single pole) in a 
fan-beam-like pattern, as illustrated in figure 12. The beam will 
point to us twice per P prec = 2.07 s, when the pole appears from 
the rim of the NS, to be called “dawn”, and disappears behind 
it, to be called “dusk”. This explains the double-peaked pulse 
profiles (figure 2, figure 7), and the ~ ±1/4 cycle offset of the 
hard X-ray pulse peaks from those in the soft X-rays (figure7b). 

As already described in section 1 and section 4.2.1, we fur¬ 
ther assume that the fan-beam emission pattern is asymmetric 
around x 3 . Then, at a certain time to of the 36 ks slip period, 
the signal will be stronger at the “dusk” pulse phase than at the 
“dawn” phase, as illustrated in figure 12 a; the pulse peak is de¬ 
layed by ~ 1/4 cycle from the average, just as in the red to 
orange modulation phases in figure 2b. At t ~ to + T/4, the 
configuration will become as in figure 12 b, so that the inten¬ 
sity will be similar between “dawn” and “dusk”; the pulse-peak 
phase will occur at about “noon” phase with zero delay, as in 
the yellow profile in figure 2b. Then, at t ~ to ± T/2, the signal 
will become strongest at the “dawn” pulse phase, as shown by 
figure 12c; the pulse peak preceds the average by ~ P pre c/4, as 
represented by the green to blue profiles in figure 2b. Finally, 
at t ~ to + 3T/4, the pulse amplitude becomes low, because the 
emission region is behind the horizon when the beam points to¬ 
wards us; this may correspond to the purple profile in figure 12 b. 
Thus, the fan-beam configuration with a broken axial symmetry 


can explain the large modulation amplitude of IE 1547.0—5408 
reaching A ~ T/4. Although the actual configuration (e.g., con¬ 
tribution from the other pole) would not be as simple as this, 
more detailed modeling of the emission pattern is beyond the 
scope of the present paper. 

How about the remaining items No.2 and No.3? In the 2009 
data of 4U 0142+61, the hard X-ray pulse profile was richer in 
fine structures, consisting of three peaks separated by a quar¬ 
ter cycle (figure If of MEA14). This is considered to be the 
reason why n = 4 was needed. Moreover, such a complex 
pulse profile of 4U 0142+61 is thought to have been strongly 
smeared out by the precession, and hence the hard X-ray pulse 
was difficult to detect without demodulation at least in the 2009 
Suzaku observation. In contrast, the pulse profiles of the present 
source are basically double-peaked. This is considered to make 
n = 2 appropriate for the present source, and make the profile 
less affected by the phase modulation. In addition, the effects 
in IE 1547.0—5408 can be approximated as a sort of pulse- 
profile change (figure 12), with the “dawn” and “dusk” inten¬ 
sities changing across the slip period, rather than a continuous 
phase shift. This property is thought to be another factor that has 
made the pulse easier to detect without demodulation. Thus, the 
differences No.2 and No.3 can be explain in a qualitative way. 

4.4 Relations to other neutron stars 

So far, a fair number of reports have been made on possible de¬ 
tections of precession from NSs, including accretion-powered 
objects (e.g., Staubert et al. 2009) and rapidly rotating NSs 
(e.g, Shabanova et al. 2001; Chukwude et al. 2003; Haberl et 
al. 2006). However, as discussed in MEA14, such a behavior in 
accretion-powered systems should be regarded as forced preces¬ 
sion, because of the strong accretion torque. The fast-rotating 
NSs can exhibit free precession, because they must be deformed 
into oblate shapes by centrifugal force. However, the precession 
would be damped on rather short time scales (Cutler 2002), be¬ 
cause an oblate body attains its minimum energy at a = 0 under 
a constant L. Therefore, these reports remain generally rather 
unconvincing. 

In contrast to these ordinary NSs, magnetars are expected 
to more ubiquitously exhibit the free precession, because their 
strong B t will deform them into prolate shapes (Ioka 2001; 
Cutler 2002; Iokai & Sasaki 2004); their centrifugal deforma¬ 
tion (oblate) is negligible, e.g., |e| < 10~ 6 , because of their 
rather long pulse periods, P = 2 — 11 s. Actually, the present 
data may favor the prolate (e > 0) case (subsection 4.2.2). In 
such a prolate magnetar, the free precession will spontaneously 
develop, because internal energy dissipation processes, with L 
conserved, will increase a on reasonable time scales (Cutler 
2002). (As mentioned in MEA14, decay in a by emission of 
gravitational waves is negligible.) Therefore, the two necessary 
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conditions, e^O and a ^ 0 , can be regarded as almost inher¬ 
ent to magnetars with strong B t . Actually, the strong pulsations 
observed from nearly all magnetars imply q^O. 

The remaining requirement for the detection of free preces¬ 
sion, i.e., the non-axi-symmetric emission pattern of the hard 
component, is considered more conditional rather than intrin¬ 
sic. One possibility is that a magnetic pole on x 3 is magneti¬ 
cally connected to a nearby local multipole, and this connecting 
region is responsible for the hard X-ray generation. Depending 
on the detailed magnetic configuration, such a radiation may 
well have the required anisotropy around, and positional off¬ 
sets from, £ 3 . Since the magnetic configuration of a magne- 
tar may change with time (Thompson & Duncan 1995), e.g., 
across active states, A can change even though T and a are 
kept constant. Actually, Kuiper et al. (2012) found significant 
pulse-profile changes in IE 1547.0—5408 across the 2009 ac¬ 
tivity episode. Similarly, the hard X-ray phase modulation of 
4U 0142+61, detected in 2009 and reconfirmed in a subsequent 
observation in 2013 (Makishima 2014), was absent in an earlier 
Suzaku observation made in 2007 (Enoto et al. 2011; see also 
MEA14). It is hence of significant interest to clarify whether the 
phase-modulation amplitude A of IE 1547.0—5408 have been 
changing along with its long-term decay after the 2009 activity. 
Future observations of this object with NuSTAR and ASTRO-H 
are very important. 


4.5 Toriodal magnetic fields 

From equation (2) and equation (17), IE 1547.0—5408 is con¬ 
cluded to be aspheric by e = Pq/T = 0.6 x 10~ 4 , which is of the 
same order as that of 4U 0142+61, 1.6 x 10 -4 (MEA14). Then, 
adopting the magnetic deformation interpretation, these NS are 
both inferred via equation (1) to have B t ~ 10 16 G. The present 
result thus gives a strong support to the suggestion by MEA14, 
that magnetars have ultra-intense B t and are hence deformed to 
a detectable level. 

While 4U 0142+61 is a rather aged magnetar with a char¬ 
acteristic age of t c = 68 kyr and Bd = 1.3 x 10 14 G, IE 
1547.0-5408 has r c = 1.4 kyr and B d = 2.2 x 10 14 G. Even 
though t c of magnetars can be systematically overestimated 
(Nakano et al. 2015), IE 1547.0—5408 is clearly a much 
younger and more active object than the other, as evidenced by 
its faster rotation and higher activity. The somewhat smaller e of 
IE 1547.0—5408 apparently contradicts to these facts, and the 
prospect that magnetars must be consuming their toroidal (or in¬ 
ternal), as well as dipole (or external), field energies. Possibly, 
the B t vs Bd ratio may scatter considerably among mangetars. 
To answer this issue, we clearly need to increase the number of 
detections of the free precession from magnetars. 


4.6 The wobbling angle a 

Finally, we consider possible values of the wobbling angle 
a. Unlike e which can be accurately determined via equation 
( 2 ), a is rather difficult to unambiguously estimate from the 
data, because its effects degenerate with those of the radiation 
anisotropy. Below, some qualitative consideration is carried out. 

On one hand, a cannot be close to 0, for the following rea¬ 
sons. First, an aligned rotator with a ~ 0 (x 3 || L) would not 
show the soft X-ray pulsation, as already noticed in section 4.4. 
Second, a would increase on a reasonable time scale (section 
4.4). Finally, the dipole field derived from the \JPP factor is 
in reality B d sin a rather than B d itself, so a ~ 0 would demand 
too high a value of B d . 

On the other hand, a cannot be close to 90°, either. This is 
because such an orthogonal rotator would exhibit a very high 
pulse fraction (unless we are observing from a pole-on direc¬ 
tion), contradicting to the fact that IE 1547.0—5408 has a par¬ 
ticularly low pulse fraction among the known magnetars. All 
what can be said at present is that a is likely to be in an in¬ 
termediate range between 0 and 90°. A more realistic estimate 
on a will be available in future when we model the hard X-ray 
emission pattern, and fit its predictions to the actual data. 

5 Conclusoin 

During the 2009 January outburst of IE 1547—5408, its 2.07 s 
pulsation exhibited a phase modulation with a period of 36 ks, in 
the 15-40 keV HXD data as well as in the XIS data at energies 
above 4 keV. The modulation is consistent with being solely 
carried by the spectral hard component, while absent in the soft 
component. These results can be understood by presuming that 
the NS in this object is exhibiting torque-free precession under 
a deformation by e = 0.6 x 10 -4 , and its hard X-ray emission 
is asymmetric around the NS’s symmetry axis. These results 
further suggest the presence of intense toroidal magnetic fields 
of B t ~ 10 16 G. 
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Appendix: Phase-randomized demodulation 

The significance of the 36 ks HXD phase modulation described 
in section 3.2.4 was evaluated through a control study using 
the same actual data themselves, rather than via a Monte-Carlo 
method. Let us explain how this was carried out. 

For this purpose, we again refer to the two-dimensional ar¬ 
ray C(j,k) (section 3.2.1), into which all the 15-40 keV HXD 
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events have been sorted. The pulse phase and the modulation 
phase are represented by k and j, respectively. The Z% calcula¬ 
tion effectively produces a folded pulse profile as 

j- i 

D{k)=Y J C{j,k ) ( 20 ) 

3=0 

and applies the Fourier transform to D(k). Here, we employ a 
rather large values of J and K, namely J = 90 and K = 360, 
to avoid binning issues. The demodulation procedure, instead, 
executes the summation by shifting the pulse phase, according 
to equation (8), as 

j 

D(k) = ^C(j,k'(j)) (21) 

3=0 

with 

k'(j) mk- K(A/P 0 h) ■ sin(27r j/J - 0) . (22) 

Let us assume, as a null hypothesis, that the pulse phase is 
intrinsically not modulated at the period T. Then, the demodu¬ 
lation procedure using equation (21) and equation (22) has two 
effects; (i) to degrade the underlying pulse coherence, and (ii) 
to add up Poisson noises of different pixels in various different 
ways. Large values of X max , like in equation (9), will appear 
when (ii) overwhelms (i). Importantly, these two effects must 
remain unchanged even when a random permeation is applied 
on the rows {j = 0,1,.., J — 1}, because all rows must be equiva¬ 
lent with respect to (i), and the Poisson noise in different pixels 
of C(j, k ) should be mutually independent when considering 
(ii). 

With the above in mind, we repeated the same demodulation 
search 1,000 times, over the same ( P,T,A,(j >) volume. In each 
trial, k'(j) in the above two equations were replaced withfc'(y'), 
where {f = 0,1,.., J — 1} is a random permutation of {j = 
0,1, ..,J — 1}. Different trials employed different realizations of 
the permutation, all with J = 90. We thus obtained 1,000 values 
of Wmax, of which the histogram is shown in figure 5a. By 
integrating this distribution from X max = +oo towards lower 
values X max , the upper probability integral has been obtained 
as in figure 5b. This gives the chance occurrence probability of 
large values of X max , under the presence of the pulsation at Po 
but the absence of any intrinsic periodicity at T. Equation 11 
provides its approximation at A' max > 35. 

This control study is considered to be rather conservative 
for the following reasons. When the phase modulation is sinu¬ 
soidal, a small difference in the initial value of <j> will not affect 
X, because At is expected to have a long coherence length in rj>. 
This makes a scan over <f> only partially independent. However, 
when the random permutation is employed, At attains a rather 
short (typically, 2n/J radian) coherence in f>, so that scans over 
4> become essentially all independent. As a result, one control 
run is expected to involve effectively a lager number of trials. 


If this effect is taken into account, the true chance probability 
would be even lower than equation (11). 
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Fig. 1. (a) Periodogams of IE 1547.0-5408 calculated using the method with n = 2, from the 2-10 keV XISO data (blue) and those of HXD-PIN in 12-70 
keV (brown), (b) The same as panel (a), but for the 15-40 keV HXD-PIN data, before (black) and after (red) applying the demodulation analysis employing the 
parameters of equation (10). In both panels, the green abscissa at the top indicates the period T which can explain a side lobe at that period in terms of its 
beat with the pulsation. The implied values of T (in units of ks) of major side lobes are also indicated in the figure. 
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Fig. 2. (a) The 15-40 keV background-inclusive pulse profiles folded at equation (4), obtained in 6 consecutive time segments of ~ 14 ks each (from red to 
purple). They are shown for two cycles, with vertical offsets, after applying running averages over 3 adjacent bins. A typical ±l<r statistical error is given at the 
top right corner. The dashed black profile at the middle represents the average, (b) The 15^10 keV folded pulse profiles accumulated over six different phases 
of the T = 36 ks modulation period, shown with offsets. Meaning of the colors is different from that in panel (a). A possible phase assignment is indicated 
by arrows, (c) The same as panel (a), but prior to the folding, the arrival time of each event was corrected with equation (8) using the parameters of equation 
(10). (d) A comparison between the 15^10 keV pulse profiles before (black) and after (red) the demodulation, both with the running averages. The former is 
the same as the black profiles in panels (a) and (b), while the latter is identical to the black one in panel (c). The background level is at about 0.27 c s _1 . 

























































































14 


Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 




Fig. 3. (a) Background-inclusive 1 5—40 keV light curves of IE 1547.0-5408 each covering the entire observation with 20 segments (j = 0,1,.., 19), derived 
in 5 different phases (red, orange, green, blue, and purple) of the pulse period of equation (4). See text for details. They are shown after applying a running 
average over the adjacent 3 bins, and their average is given in solid black line. Typical ±io error bars are given at the bottom right corner, (b) The same as 
panel (a), but sorted into 5 phases of a dummy period of 2.0 s. (c) Power spectra of the five light curves in panels (a) and (b), in red and blue lines, respectively. 
Both are normalized to the power expected from the Poisson fluctuations. The black dashed line represents an ensemble average over 20 control studies, 
each similar to panel (b). 
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Fig. 4. Results of the demodulation analysis on the 15—40 keV HXD-PIN data, assuming a periodic phase modulation described by equation (8). (a) A two- 
dimensional color map of X = z\ (maximum as P is varied) on the ( A , 4>) plane, assuming T = 36.0 ks. (b) Projection of panel (a) onto the <f> axis, (c) 
Projection onto the A axis, (d) T-dependence of the z\ peak value, 7f pk (T), presented in black. Green shows a reference result using the 300-307 keV 
HXD-GSO data. The corresponding Fourier wave number, for the total time span of 86 ks, is shown in blue at the top. 
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Fig. 5. A histogram of the values of X max (panel a), and the associated upper probability integral (panel b), found in a series of 1,000 simulation trials each 
equivalent to figure 4 but randomizing the 36 ks modulation phase. As described in Appendix, it was derived from the data themselves. In both panels, a black 
arrow indicates equation (9). 
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Fig. 7. Effects of the demodulation on the 10-14 keV XISO data. Equation (14) and T — 36.0 ks are employed, (a) Periodgrams, before (black) and after (red) 
the demodulation. The plot is more expanded in the period range than those presented so far, to show finer details around the pulsation, (b) The 10-14 keV 
pulse profiles folded at P = P 0 with the colors corresponding to those of panel (a). Running averages over the adjacent 3 bins are applied. The phase origin 
is the same as in figure 2. Green (with right ordinate) shows the 1-10 keV profile. 
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Fig. 8. (a) The same A-dependence of A' as figure 4c and figure 6c, but for the 7-10 keV XISO data. T = 36 ks is fixed, while 4> is left free, (b) In 4-7 keV. 
(c) In 1^1 keV. (d) The same as figure 2b but for the 7-10 keV XISO data. 
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Fig. 9. The 36 ks phase modulation amplitude A, obtained at various energies with the HXD and the XIS. The modulation phase is not specified. The dashed 
curve represents equation (16). 
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Fig. 10. The 15-40 keV pulse profiles in the 1st (red) and 2nd (blue) 36 
ks, both demodulated using the parameters in equation (10) and folded at 
equation (4). Their average, shown in dotted black lines, can be slightly 
deferent from the red data in figure 2d because of the exclusion of some 
fraction of the data (see text). 
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Fig. 11. All the 15—40 keV periodograms obtained in computing figure 4, 
superposed onto a single plot. The two periodograms in figure 1b are both 
included here. 
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Fig. 12. A possible free precession scenario of an oblate NS, as seen from above onto its angular momentum vector L whose tip is represented by a black 
star. The observer is viewing from the left. The NS’s symmetry pole is represented by a small blue circle, and its precession locus by a larger dashed blue 
circle with a blue arrow. Green ellipses show hard X-ray emission pattern, which is assumed to be fan-beam like. The panels (a), (b), (c), and (d) represent 
slip-period epochs of to, t 0 + T/4, t 0 + T/ 2, and to + 3T/4, respectively. A and B indicate “dawn” and“dusk”, respectively. Lower panels show the expected 
hard X-ray pulse profiles. 







